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Future metrology standards will be partly based on physical quantities computed from first princi¬ 
ples rather than measured. In particular, a new pressure standard can be established if the dynamic 
polarizability of helium can be determined from theory with an uncertainty smaller than 0.2 ppm. 

We present calculations of the frequency-dependent part of this quantity including relativistic effects 
with full account of leading nuclear recoil terms and using highly optimized explicitly correlated ba¬ 
sis sets. A particular emphasis is put on uncertainty estimates. At the He-Ne laser wavelength of 
632.9908 nm, the computed polarizability value of 1.391 811 41 a.u. has uncertainty of 0.1 ppm that 
is two orders of magnitude smaller than those of the most accurate polarizability measurements. 

We also obtained an accurate expansion of the helium refractive index in powers of density. 


Some physical quantities, for example, properties of 
the helium atom and interaction energies of helium 
atoms, can now be computed from first principles with 
precision rivaling and sometimes exceeding the best ex¬ 
perimental determinations [1-3]. Therefore, quantities 
of this type can be used in establishing metrology stan¬ 
dards. One example is a possible standard of temper¬ 
ature based on acoustic gas thermometry [4]. Another 
example is a pressure standard based on optical interfer¬ 
ometry [5]. The current pressure standard dating back 
more than 300 years is realized by mercury manometers 
and can not be further improved. Also, the reference 
manometers are far from portable: 3 m high and con¬ 
taining 250 kg of mercury, a substance banned due to 
its toxicity. Since pressure is one of the most widely 
measured properties, in applications ranging from man¬ 
ufacturing of semiconductor chips to air-traffic control, 
a new pressure standard would significantly impact both 
technology and everyday life. The proposed standard [5] 
obtains pressure from the formula [6] 

_n^-l ikT 
^ n? -\- 2 47r(a + x) ’ 

where n denotes the index of refraction of helium gas, k 
the Boltzmann constant, T the temperature, a the dipole 
polarizability, and x the diamagnetic susceptibility of he¬ 
lium. To account for nonideality of helium gas, one has to 
include some small terms on the right-hand side depend¬ 
ing on dielectric and density virial coefficients [6]. The 
essential part of the new standard is the determination 
of n with an uncertainty of 0.2 ppm via interferometric 
measurements of a variable-length cavity filled with he¬ 
lium and comparing to measurements in vacuum. The 
product kT, currently known with an uncertainty of 0.9 
ppm [7] near the temperature of the triple point of water, 
is the subject of active research and a reduction of this 
uncertainty can be expected in near future. Since x is five 
orders of magnitude smaller than a, it can be computed 


using the non-relativistic wave function from the expres¬ 
sion X = —e^{r^)/3TOeC^, where e and We are the electron 
charge and mass, c is the speed of light, and (r^) is the 
average square of the electron-nucleus distance. Also the 
virial coefficients are known accurately enough from the¬ 
ory [3, 8]. However, a cannot currently be measured with 
uncertainty lower than 0.2 ppm, so the standard clings 
upon theory being able to achieve such accuracy. This 
Letter describes calculations of a from first principles. 
We also discuss how n obtained purely from theory can 
be used to calibrate refractometers or to correct errors in 
interferometric length measurements [9]. 

Since the radiation frequency of interest, 632.9908 nm 
[10], is much smaller than the lowest resonance, the fre¬ 
quency dependence of a can be efficiently calculated 
from the power series expansion, a(oj) = ag -h a 2 + 
04 U!^ -h • • •, where ag is the static dipole polarizabil¬ 
ity. The coefficients Ok, k > 0, describing the frequency 
dependence of polarizability, will be referred to as the 
(polarizability) dispersion coefficients. We shall use the 
atomic units throughout (we never use reduced atomic 
units), in particular, Og, where ag is the bohr radius 
k ?as the unit of polarizability and the inverse 
of the atomic unit of time, tg = /me,e^, as the unit of 

frequency. For light systems like helium, each Ok can be 
expanded in powers of the fine structure constant 1/c, 
where c=137.0359991 [7] is the speed of light expressed 
in atomic units, au = -\- + • • •, be¬ 

ing proportional to 1/cK We shall refer to as the 
relativistic contributions. The terms sfc., are 

due to radiative as well as higher-order relativistic effects 
predicted by quantum electrodynamics (QED). 

The nuclear mass dependence of the nonrelativistic po¬ 
larizability an \ can be taken into account exactly, but 
for the relativistic and QED contributions one has to 
use an expansion in powers of the ratio of me to the 
nuclear mass rria, i.e., in powers of 1/M = me/ma = 




2 


1/7294.2995361. Since 1/M is of the order of 10 
keeping the linear term is entirely sufficient and such 
contributions can be represented in the form = 
I > 2, where are computed with the in¬ 
finite nuclear mass and are corrections of the order 
of l/(Mc*), referred to as the recoil corrections. These 
recoil corrections are expected to be negligible except for 
the static ones and and, possibly, for 

For comparisons with experiments, it is convenient to 
convert frequency to wave length A = 27rce^/fia; 

a(A) = Ao+A 2 A“2 -b ^4 A"^ H-. (2) 

When a(A) remains in atomic units and A is measured 
in nm, the relation is A^ = f^o-k, f = 27rcao/nm = 
45.56335253 (with ao=0.05291772109 nm ). 

For ao contributions, our more accurate values are con¬ 
sistent with Ref. 1 to all digits published except for the 
term describing the electric-field dependence of Bethe’s 
logarithm and for . We computed the former term 
using a new method since Ref. 1 was the only source for 
this quantity. The contributions of this component from 
the two calculations differ only marginally, by 0.011 ^Uq. 
The term was estimated in Ref. 1 by the contri¬ 
bution from the simple one-loop expression [11] and the 
uncertainty of this term was assumed to be 40%. Later, 
it was shown in Ref. 12 that the error of one-loop ap¬ 
proximation applied to the excitation energies of helium 
is only about 5%. Therefore, we reduced our estimate 
from 40% to 25% or 0.14 /toq, which we believe is still 
conservative. 

The dispersion coefficients ak {k = 2, 4, 6) were cal¬ 
culated thus far only by Bhatia and Drachman (BD) 
[13, 14]. However, these authors did not provide any 
estimates of the uncertainties. Their relativistic con¬ 
tributions do depend on the nuclear mass but the re- 
( 21 ) 

coil effect, a/ , was not correctly taken into account, 
vide infra. Furthermore, the Ak coefficients were incor¬ 
rectly converted from the reduced Rydberg units: the 
factor (1 -b me/ma)^, appearing in the correct conver¬ 
sion formula, was erroneously replaced by its square 
{l+me/rriaY^. 

At the nonrelativistic level of theory, a{oj) of an atom 
in a quantum state ip is defined by the standard polar¬ 
ization propagator expression 

a{uj) = {'tp\zTZ{oj)z\\p) -b {ii\zTZ{-uj)z\'il)), (3) 


where z = Zi+ Z 2 , with Zi denoting electron coordinates, 
and 'R,{ijj) = Q {QH — E + uj)~^ is the resolvent of the 
atomic Hamiltonian H, with Q = 1 — P = 1 — \i/){-tpl and 
E being the energy of state ip. For the helium atom 


H = 




2M 


(Vi+V2f 


2 
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2 
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1 

ri2 


( 4 ) 


TZ{uj) satisfies the identity = TZ — loTZTZ{lo), where 
TZ = Q {QH — E)~^ is the static (reduced) resolvent of 


H. Iterating this expression and inserting it into Eq. (3), 
one obtains 

af =2{iP\zTZ’^+^z\'tP). (5) 

To account for the leading relativistic contributions of 
the order of 1/c^ assuming infinite nuclear mass, we add 
to the Hamiltonian of Eq. (4) the perturbation from the 
Breit-Pauli Hamiltonian [15] obtaining 

= “4 (V’o|Bi77oz77.ozV’o) - 2 {'4}o\zTZqBiTZqziI}q), 


a^2°'^ = -4 (V'o|Bl7^oz7^oZ■^/’o) - 4 {\po\zTZoBiTZlz'tj)o) 
-2 {iPq\zTZIBiTZIz%Pq), 

af = -4 - 4 {MzT^oBiTZlz'i/Q 

-4 ('ipojzHo^i'H^zipo} - 2 {'4 }q\zTZIBiTZIziI}o), 


= -4 {'ipo\BiTZozTZlz%pQ - 4 
-4 {iPo\zTZIBiTZIz4)q) - 4 
-2{iPo\zTZtBintz4,o), 

where 

= -^(V^ + V^) + ^[5{rQ + S{r2)] + ^S{r,2) 

+ ^[Virf2V2 -b (Viri2)rf2^(ri2V2)], (6) 

Bi = Bi — {'ipo\Bi'ipo), and the quantities with subscript 0 
are analogous to those defined above but for infinite-mass 
nonrelativistic Hamiltonian. 

( 21 ) ( 21 ) 

The relativistic recoil term ag ^ is equal q;q (B 2 ) + 
a''g^\HiBi), where Hi = -^{^i + ^ 2 ^ and 


B. = 


Men 


Vir;^ ^Vi -b (Viri)ri ^(»’iVi) -b Vir^ V 2 


+ (Viri)ri ^(riV2) -b V2r2 ^Vi + (V2r2)r2 ^(r2Vi) 
-b V2r2^^V2-b (V2r2)r^^(r2V2)j. (7) 


The two components are given by 

a''o^\B2) = -4 {tpo\B2TZozKoztpo)-2 {tpo\zTZoB2TZozipo) 

( 8 ) 

where B 2 is analogous to Bi and 
al^^\HiBi) = 

4 [{ipolzTZoZTZoHiTZgBilpg) + {/^glzTZQzTZQBiTZgHilpg) 
-b (iiolzTZoHiTZozTZoBi-tpo) + (tpolzTZoBiTZozTZoHi/jo) 
-b {i^olz'JZoHiTZoBiTZoZ'tjjo) + (f/'oj-ffiTT.o^T^o^TT.o.Bi'i/o) 

- {'ipo\z'IZoZ'ipo){'ipo\HiTZlBi'tpo) 

-{'tpo\zTZlz'ipo){'ipo\HiTZoBi'tpo)] , (9) 

where Hi = Hi - (ipolHitpo). 
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The correction is very small and can be computed 
using a finite difference expression 

^ B 2 ) + a^2°\Ho ^ H) - (10) 

valid to the order of l/(M^c^), where Bi —> B 2 means 
that operator Bi in the expression for should be 

replaced by B 2 and similarly Hq —>■ H means that quan¬ 
tities computed with the Hamiltonian Hq should be re¬ 
placed by those computed with H. 

To evaluate accurate representations of the he¬ 
lium ground-state wave functions ipo and ip were obtained 
by minimizing the conventional Rayleigh-Ritz functional 
for the Hamiltonians Hq and H, respectively. The auxil¬ 
iary functions were obtained recursively from Hylleraas- 
type functionals 

= {^\Ho -Eo + Po|0) - (11) 

for = TZqZiPq and 

= {mQ-EQ + PQ\^)-2{^\(l-PQ)z4l'Q^), ( 12 ) 

for = TZqz4>^q\ and analogous functionals obtained 

by dropping all the subscripts 0 for ())(”) = HP zip and 
^(") = TZz(p^'^\ The trial functions used in all minimiza¬ 
tion processes were expanded in bases of Slater geminals 

N 

0 = (1 + V 12 ) r (n, ra) , (13) 

where 7^12 is the transposition operator whereas 
y(ri,r 2 ) = Zi in calculations of and and 
y(ri,r 2 ) = 1 otherwise. One may note that the func¬ 
tions and ip^^^ contain also a D component, but it 
does not contribute to matrix elements that are needed. 
The linear coefficients were obtained by solving the ap¬ 
propriate set of linear equations, while to determine the 
nonlinear parameters we employed two strategies: the 
full optimization (FO) and the stochastic optimization 
(SO). In the latter case, the parameters a^. Pi, 7 * are 
pseudo-randomly generated from a box with optimized 
dimensions. We used two boxes to model the short-range 
and medium-range asymptotics of the wave functions. To 
eliminate possibilities of numerical errors, the FO and 
SO based codes (including the integral and linear alge¬ 
bra routines) were programmed entirely independently 
by different members of our team. 

The contributions for k = 0, 2, 4, 6 were com¬ 
puted for several values of N, up to 600 (800) in the 
FO (SO) approach, both optimizations giving at least 11 
convergent digits, with FO converging faster. Our results 
agree to 9, 8 , 4, and 7 digits, respectively, with the val¬ 
ues obtained by BD [13]. Using the SO procedure, we 
also calculated: 4°^ =4.39500532(1), a® =6.7725956(1), 
42^=10.622083(1), and a^ 4 ^= 16 . 86118 (l) fxal 


For the relativistic contributions fc = 0, 2, 4, 6 , 

the convergence is much slower than in the nonrelativistic 
case. This is due to the fact that we use nonrelativistic 
functionals which are sensitive to wave function values in 
different regions of the configuration space than the rela¬ 
tivistic operators (these operators are too singular to be 
used in optimizations). The SO procedure leads now to a 
faster convergence than FO since randomly chosen expo¬ 
nents cover the space more uniformly than FO exponents. 
Thus, we used the SO results as our recommended values 
and in estimates of uncertainty. Nevertheless, the agree¬ 
ment to 6 , 5, 3, and 3 digits, respectively, between the 
two sets of results is more than sufficient for the present 
purposes. Our values are substantially more accurate 
than those of BD [14], with agreement to only 2, 2, and 
1 digit, respectively (BD did not compute For 

a.^Q^\ our results are consistent with, but significantly 
more accurate than calculations of Refs. [1, 16, 17]. The 
relativistic recoil contribution ag^^^=—0.0935(1) /iOg. Its 
smallness results from some cancellation of its compo¬ 
nents, ag^^^(i7ii?i) and a'Q^\B 2 ), equal to 0.1559 and 
—0.2494 ^Oq, respectively. The contribution is 

equal to —0.144(1) ^Og, so it is virtually negligible. 

It should be pointed out that the relativistic contri¬ 
butions computed by BD [14] depend on the nuclear 
mass and, strictly speaking, should not be compared 
with our, nuclear-mass-independent contributions 
This is because these authors incorrectly assumed that 
the individual terms in the Breit-Pauli Hamiltonian are 
proportional to (inverse) powers of the reduced electron 
mass rather than the real mass. Therefore, although the 
nuclear-mass-dependent part of their relativistic contri¬ 
butions is of the order of l/(Mc^), it differs from the 
ap (BiHi) part of the true recoil correction. Addition- 
ally, BD completely neglected the contribution ap ' (52 )■ 
Thus, their relativistic contributions cannot be viewed as 
approximations to Since the effects of the 

order of 1/{M(P) are very small, the differences between 
our relativistic contributions and those of BD are mainly 
due to the differences in basis sets used in the calcula¬ 
tions rather than to the treatments of the nuclear mass 
dependence. 

After correcting the units conversion error in Ref. 14 
as discussed earlier, the coefficients computed by BD 
agree with our values to 5, 6 , 4, and 5 digits for k = 
0, 2, 4, 6 , respectively. For fc = 0, the discrepancy is 
mainly due to the 1/c^ terms not considered by BD. The 
reasons for the relatively low accuracy of A 4 are unclear. 
Due to the smallness of the relativistic contributions to 
Ak, the overall agreement is good despite the fact that the 
relativistic contributions from BD work are significantly 
less accurate than ours. 

In Table I, we present the dynamic polarizability of 
^He. In addition to the contributions discussed earlier. 
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TABLE I. Dynamic polarizability of ^He [aq] at A = 632.9908 
nm. 


static 

nonrelativistic 

l/c^ “ 

l/c^ 

llc^ 

finite nuclear size“ 
total 

1.383 809 98641(1) 
-0.000 080 4534(1) 
0.000 030 655(1) 
0.000 000 56(14) 
0.000 000 0217(1) 
1.383 760 77(14) 

A-’" 

nonrelativistic 

0.007 995 7979(1) 


relativistic'* 

-0.000 0001721(1) 


total 

0.007 995 6258(1) 

A''* 

nonrelativistic 

0.000 054 8363(1) 


relativistic 

0.000 000 00014(1) 


total 

0.000 054 8364(1) 

A-'’ 

nonrelativistic 

0.000 000 4076(1) 


relativistic 

0.000 000 0000(1) 


total 

0.000 000 4076(1) 

A-“ 

nonrelativistic 

0.000 000 0032(1) 

a{\) - 

a(0) present“ 

0.008 050 8730(1) 


BD^ 

0.008 050 871 

total 

present 

1.39181164(14) 


BD» 

1.391 780 800 


“Includes the recoil correction of the order of l/(Mc^) equal 
to —0.000 000 0935(1). *’From Ref. 1 except for the 
contribution from the electric field derivative of the Bethe 
logarithm equal to 0.000 000 182(1) [18]. “ Computed 
adding the correction term (4/3)7rr„ [(5(^1) + 5(r2)] to H, 
where Va = 1.676 fm is the nuclear radius. ‘^Including the 
recoil correction of the order of lf{Mc^) equal to 
—0.000 000 00075(1). “The contribution of the A“^° term, 
amounting to 2.5x10“^^, is negligible. -^Calculated using 
correctly converted Ak constants. Equation (15) of Ref. 14 
gives 0.008 052 951, i.e., 0.03% error resulting in 1.7 ppm 
error in the total value of a(632.9908). ®Using the static 
value of BD equal to 1.383 729 929. 

we included the effect of finite nuclear size which is almost 
negligible. The dispersion part of a(632.9908), i.e., the 
contribution explicitly dependent on wavelength, agrees 
to 6 significant digits with the result of BD (after conver¬ 
sion errors are corrected) due to the high, eight-digit ac¬ 
curacy of BD’s contribution. However, the total po¬ 
larizability obtained by us differs significantly, agreement 
to 5 digits and discrepancy of about 22 ppm, from BD’s 
result. As already discussed, this difference is mainly 
due to the QED effects neglected by these authors. The 
second source of the difference is our significantly im¬ 
proved value of the static relativistic component. The 
uncertainty of our recommended value of a(632.9908) 
amounts 0.14 /rog, i.e., about 0.1 ppm. This accuracy 
is sufficient for the purpose of the new pressure standard 
but one should ask if any neglected effects could con¬ 
tribute above the uncertainty estimate. The potential 
candidates are the QED recoil correction ag of the or¬ 
der of l/(Mc^), the QED contribution to the polarizabil¬ 
ity dispersion of the order of 1/c^ and, finally, and 
probably most importantly, the remaining, other than 
one-loop contributions to of the order of 1/c^. We 


believe that such neglected contributions should not con¬ 
tribute more than 0.1 ppm but are investigating such 
terms. 


TABLE 11. Virial expansion of refractive index. Ur = ^a„, 
bn, and Cn are in units of cm®/mol, cm®/mol^, and cm®/moE, 
respectively, and A is in nm. 1 cm®/mol = 11.205 8721 Aq. 


ar(0) 

present 

0.517 246 21(6)“’'’ 


exp. [6] 

0.517 245 5(47) 

a^ (632.9908) 

present 

0.520 255 64(6)“’“* 


exp. [19, 20] 

0.5213(1) 


exp. [21] 

0.522 0(3) 

a^ (546.2268) 

present 

0.521297 25(6) 


exp. [22] 

0.52157(15)“ 

6n(0) 

present, 273.16 K 

0.0245(2)-'' 

6„ (632.9908) 

present, 273.16 K 

0.0238(2)-*' 


present, 302 K 

0.0184(2/ 


present, 323 K 

0.0151(2/ 


exp. [20] 

0.000(15/ 

c„(0) 

present, 273.16 K 

-0.93(25)-*' 


“ a„ = 8.694 292 2(9) ag or 0.775 869 31(8)(4) cm^mol, 
where the second uncertainty originates from the Avogadro 
constant 6.022 141 29(27) x 10^®. ** Computed using y = 
—0.000 021 194(1) ag [23], the uncertainty reflects the 
estimated size of relativistic contributions. “ a„ = 

8.744 8758(9) ag = 0.780 383 35(8)(4) cm^mol. We 
neglected the frequency dependence of x, a relativistic effect 
expected to be very small. “Inferred from measured value of 
n-1 = 34.895(10)xlO”"' at r=273.16 K and p=101.325 
kPa. Computed using be=—0.0978(2) cm®/mol [24] 
(uncertainty estimated based on comparison with Ref. 8) 
and Ce = — 1.34(36) cm®/moE [6, 25], r=273.16 K. For other 
T: be(303) = -0.1065(2) cm^mol and &e(323) = 

—0.1107(2) cm®/mol [24[. ® Computed from a^ = 0.5213(1) 
cm®/mol and br = |6n — ^cir = —0.068(10) cm®/moff 
measured in Ref. 20, same value obtained for both 
temperatures. 

The virial expansion for the refractive index can be 
written as 

n = 1 + GnP + bnP^ -|- CnP^ + • • • , (14) 

a„ = 27r(a-h x), (15) 

bn = 2Tr{abs + + x^m + + 27rQ;x), (16) 

c„ = 27r(Q;C£ -h b^ + (17) 

where p is density, bg and Cg are the dielectric virial co¬ 
efficients and b^ is the magnetic permeability virial coef¬ 
ficient. We have written down the term xbfi in Eq. (16), 
but we will neglect it in numerical calculations since x 
is about five orders of magnitude smaller than a and b^ 
(unknown) is expected to be at the most of the same or¬ 
der as bg. We completely ignored the magnetic part of c„ 
in Eq. (17). After Eq. (14) is squared, it becomes con¬ 
sistent with Eq. (4) of Ref. 6 within the terms included 
there except that the factor of 2 is missing in front of the 
Ag b p^ term. Equation (14) can be easily solved for p and 
the resulting formula can be used for a determination of 
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density, or, when combined with the virial equation of 
state, also for a determination of pressure. 

The virial coefficients are presented in Table II. The 
agreement with the measurement of Schmidt et al. [6] is 
excellent, to within I.4±9.1 ppm. Note that the authors 
of Ref. 6 reported the value of Ur with a subtracted mag¬ 
netic contribution dTry/S = —0.0000080 cm^/mol, which 
was added back in Table II. The agreement with mea¬ 
surements at 632.9908 nm [19-21] is, however, poor, as 
noticed earlier comparing with older theoretical results 
by BD [14] and by Stone and Stejskal [9]. The disagree¬ 
ment with the measurement of Leonard [22] is smaller, 
only about twice the experimental uncertainty. The ap¬ 
parent better agreement of theory with this experiment 
(within la) found in Ref. 14 was due to the neglect [14] 
of the nonlinear dependence of density on pressure. The 
values of and c„ presented in Table II have uncer¬ 
tainties due entirely to uncertainties of bg and Cg. The 
third and fourth term in Eq. (16) make negligible contri¬ 
butions and there is substantial cancellation between the 
first two terms. The experimental determined from 
the values measured in Ref. [20] is consistent with zero, 
which is almost within the combined uncertainties. From 
our data and from the bg{T) data of Ref. 24 we predict 
that bn will vanish only around 415 K. Since is small 
at T=273.16 K and higher temperatures, its accuracy is 
sufficient to predict n -1 with a 1 ppm uncertainty for 
pressures up to 10 MPa. 
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